Free Random Levy Matrices 
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Using the theory of free random variables (FRV) and the 
Coulomb gas analogy, we construct stable random matrix en- 
sembles that are random matrix generalizations of the classi- 
cal one- dimensional stable Levy distributions. We show that 
the resolvents for the corresponding matrices obey transcen- 
dental equations in the large size limit. We solve these equa- 
tions in a number of cases, and show that the eigenvalue dis- 
tributions exhibit Levy tails. For the analytically known Levy 
measures we explicitly construct the density of states using 
the method of orthogonal polynomials. We show that the 
Levy tail-distributions are characterized by a novel form of 
microscopic universality. 



I. INTRODUCTION 

There is a wide and growing interest in stochastic pro- 
cesses with long tails in relation to self-similar phenom- 
ena. Long tails cause the 2D random walk to be at- 
tracted to Levy stable fixed points with ubiquitous phys- 
ical properties such as intermittent behavior and anoma- 
lous diffusion. Physical examples include charge carrier 
transport in amorphous semiconductors, vortex motion 
in high temperature superconductors, moving interfaces 
in porous media, spin glasses and anomalous heat flow 
in heavy ion collisions. Levy distributions have also find 
applications in biophysics, health-physics and finances. 
Their importance stems from their stability under con- 
volution, i.e. the sum of two Levy distributed random 
variables follows also a Levy distribution. Levy stability 
for scale-free processes is the analogue of Gauss-stability 
for scale-dependent processes with finite variance. 

In the realm of complex and/or disordered systems, 
the theory of random matrices plays an important role in 
differentiating noise from information. Also, it allows for 
a generic analysis of complex phenomena in the chaotic 
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regime using random matrix universality. So far, there 
does not seem to be a developed theory of random Levy 
matrices, with the exception of j^] which is not based on 
FRV. The reasons are two-fold: Levy distributions are 
usually defined by their characteristic functions, mak- 
ing their probability density functions (pdf) usually user- 
unfriendly. Second, Levy distributions do not have finite 
second and higher moments, making standard techniques 
of random matrix theory usually of little use. 

The aim of this paper is to provide an appropriate gen- 
eralization of the classical stable one-dimensional Levy 
distributions to the random matrix setting. We use re- 
sults from the theory of free random variables, to formu- 
late a theory of stable Levy random matrices. In sec- 
tion 2, we use the explicit results from Bercovici and 
Voiculescu to show that the resolvent for free random 
variables obeys certain transcendental equations. In sec- 
tion 3, we argue that the formal result of FRV could be 
explicitly constructed using large matrices. We use the 
Coulomb gas analogy to construct the Levy measures. 
In section 4, we show how the orthogonal polynomial 
method can be applied to Levy random matrices, and 
demonstrate a novel feature of universal scaling in the 
tails. Our conclusions are given in section 5. 



II. RESOLVENTS OF FREE LEVY VARIABLES 

The fundamental problem in random matrix theory is 
to find the distribution of eigenvalues \i in the large N 
(size of M) limit, i.e 



P(A) 



\ 1 

^5(A-A0) = - (Tr<5(A-M)) 
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\i=l 



where the averaging is carried using some pertinent mea- 
sure 



-JVTry(M) 
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The potential V{M) does not have to be a polynomial 
or to have even a power series expansion, as usually as- 
sumed. Indeed, we will show below that for Levy distri- 
butions the potential is not analytic in M . In general, it 
is convenient to introduce the Green's function 
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z~ M 



(3) 



The eigenvalue distribution follows from the discontinu- 
ity of G{z) along the real axis, i.e. p{X) = — ImG(A + 
iO) /tt. In this paper we are interested in stable random 
matrix ensembles i.e. ensembles in which the eigenvalue 
distribution of the sum M = Mi +M2, is the same (up to 
a shift and/or rescaling) as for each individual matrix Mi 
and M2. If the second moment exists, then the ensemble 
is Gaussian. In the opposite case there exist, however, 
alternate ensembles which are the matrix analogues of 
the classical one-dimensional Levy distributions. 

The notion of addition laws of the type M = Mi + M2 
has been analyzed in the context of the theory of free ran- 
dom variables (FRV), a generalization by Voiculescu H] 
of classical probability theory to a noncommutative set- 
ting. A novel twist in the theory of FRV came from the 
realization, that the abstract concepts of operator alge- 
bras and free random variables can have an explicit re- 
alization in terms of large random matrices. Hence FRV 
techniques provide a novel and powerful way of analyzing 
the spectra of random matrices In this framework 

the determination of stable eigenvalue distributions can 
be performed algebraically in a very general setting. This 
has been done by Bercovici and Voiculescu ||^. However 
the drawback of such an abstract approach is that the 
explicit random matrix ensembles corresponding to these 
distributions are unknown. The aim of this paper is to fill 
this gap, as the knowledge of explicit stable random ma- 
trix ensembles with power law tails might be interesting 
for the various applications we have cited. 

The remarkable achievement by Bercovici and Voi- 
culescu 1^ is an explicit derivation of all R-transforms 
(defined by the equation R{G{z)) = z — 1/G{z) where 
G{z) = (1/(2: — A))) for all free stable distributions, with- 
out recourse to a matrix realization. Indeed, Bercovici 
and Voiculescu have found that R{z) can have the trivial 
form R{z) = a or : 



R{z) 



bz""- 



(4) 



where 0<a<2,aisa real shift parameter, and 6 is a 
parameter which can be related to the slope a, skewness 
/3, and range 7 of the standard parameterization of stable 
distribution I0,pi : 



7 e*(f-i)(i+^)'^ for 1< a < 2 



7 e 



T+f (1+/3H fQj. < a < 1 



In the marginal case: a = 1, R{z) reads : 
R{z) = a - ^7(1 + f])-^\njz 

TT 



(5) 



The branch cut structure of R{z) is chosen in such a way 
that the upper complex half plane is mapped to itself. 
Recalling that R — z — l/Gin the large N limit, one finds 
that for the trivial case R{z) = a, the resolvent : G{z) = 



{z — a)~^ and the spectral distribution p(A) = 5{X — a). 
Otherwise, on the upper half-plane, the resolvent fulfills 
an algebraic equation 

bG^'iz)- {z- a)G{z) + l^Q , (6) 

or in the marginal case (a = 1) : 

z-a + i-i{l+l3)^G{z) + "^Giz) \n-fG{z) - 1 = 0. (7) 

On the lower half-plane G{z) = G{z) The solution of 
the latter equation will not be discussed here, except for 
/3 = for which it simplifies to 



G{z)^ 
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z — (a — 17) 



(8) 



Thus in this case the spectral density has the form of a 
Cauchy distribution : (Levy with a = I) : 



TT (A — a)2 + 7^ 



(9) 



For a = 1/4, 1/3,1/2, 2/3, 3/4, 4/3, 3/2 and 2 the alge- 
braic equation m) are exactly solvable. The correspond- 
ing spectral functions are given by appropriate solutions 
of the quadratic, cubic ( Cardano) or quartic (Ferrari) 
equations. The case a = 2 corresponds to the known 
Gaussian case. We note that the spectral density for 
slope parameter 1/a follows from 



(10) 



for 1 < a < 2 and /3' = (a - 1) - (2 - a)P §, again 
by analogy to the well known duality relations for one- 
dimensional stable distributions [Q. 

Despite the fact that one-dimensional Levy distribu- 
tions are unknown explicitly (the exception being the 
three cases a ~ 1/2,1,2), the corresponding resolvents 
for free random Levy matrices obey simple multivalued 
equations. Eq. ^) has an asymptotic solution for /3 = 0, 
7=1 and A — + 0, 



p(A)^i(l-^ 



and one as A 00, 
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(11) 



(12) 



The scaling connecting the distribution with unit scale 
to the one with scale 7, 



PyW =7 " p(7"A) 



(13) 
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III. RANDOM MATRIX STABLE LEVY 
ENSEMBLES 

In our quest for a random matrix realization of the 
above free stable distributions we will look for hermi- 
tian 1^ random matrix ensembles with the measure (|^) , 
such that the corresponding p(A) coincides with the large 
N limit of the mean eigenvalue distribution. From the 
theory of FRV such ensembles will be automatically sta- 
ble with respect to addition. Here the potential V{M) 
is not restricted to be a polynomial (or a power series) 
but can involve logarithms etc. Note that in particular 
we require unitary (orthogonal) symmetry M — > UMW 
of the ensemble. 

The standard procedure of diagonalizing M UAW 
and integrating out U gives rise to the standard joint 
probability distribution for the eigenvalues 



i<j 



(14) 



For large N, the corresponding partition function can 
be analyzed in terms of the Coulomb gas action with a 
continuous eigenvalue distribution p(x) as originally sug- 
gested by Dyson M. Specifically, 



Sjp) 

iV2 



dXpiX)ViX) - / dXdX In |A-A |p(A)p(A ) , 



Two functional differentiations yield, 



V'{X) = 2PVy" 



dX 



P(A) 
A - A' 



(15) 



(16) 



where PV denotes the principal value of the integral. The 
knowledge of ^(A) plus the boundary conditions on the 
singular integral equation allows to find the spectral func- 
tion [|j . Conversely, from the knowledge of the spectral 
function we can deduce the shape of the potential, hence 
the weight of the probability distribution for the matrix 
ensemble. This is the route for implementing the FRV 
results on random matrix ensembles as we now show. 

Using the formula l/(A±ie) = PV j^i7TS{X) we get a 
formula for V^'(A) in terms of the real part of the resolvent 
on the cut V'{X) = 2 Re G(A). We can now use the resol- 
vents found in section 2 to reconstruct the potential and 
hence to explicitly define a random matrix realization of 
the free stable Levy ensembles. Below, we reconstruct 
the pertinent measures for the Cauchy (a — 1,(3 — 0) 



'^The other two possibilities, symmetric and skew-symmetric 
ones, could be considered along the lines presented in this 
work. 



and Smirnov {a = 1/2, /3 = ±1) distributions. In the 
cases a ^ 1/4,1/3,1/2,2/3,3/4,4/3,3/2 such construc- 
tion can be performed numerically. General results can 
be achieved without these specifics as we will show below. 

In the case of the Cauchy distribution, elementary inte- 
gration gives V^'(A) = 72A/(A2 + 6^), hence the potential 
V{X) — 7ln(A^ -I- 5^). In the case of the Smirnov distri- 
bution, the Green's function follows from 



—I 

7^ 



1 

G 



which can be easily solved to yield 



G(z) = 



2z - 1 - i^/Az^ 
2^2 



(17) 



(18) 



Evaluating the real part gives ^'(A) = 2/ A — 1/A^, so 
that 



V{M) = exp ( -iV Tr ( ^ 2 In M 



(19) 



The spectral function follows from (3) and (21) and reads 
p{X) = l/(27r)\/4A- 1/A2. 

In general, the asymptotic form of the potential for 
large eigenvalues reads 



V{X) = 21nA 



2 - Re 6 — 

a A" 



(20) 



In all Levy cases the 2 In A contribution in the potential 
is fixed, and generate in the measure a determinant with 
also fixed power —2N. A deviation from 2 can be shown 
to lead to a finite support of eigenvalues. The coefficient 
of the second term in the potential can vanish in some 
notable cases like e.g. for the Smirnov ensemble. In the 
next section we will analyze in greater detail the stable 
Levy random matrix ensembles defined above. 



IV. ORTHOGONAL POLYNOMIALS 

Both the Cauchy and Smirnov ensembles can be stud- 
ied analytically at finite TV using the orthogonal poly- 
nomial method. It turns out that the above ensembles 
have some unexpected new features which do not appear 
in the classical case when V{M) is a polynomial in AI. 
The orthogonal polynomials associated with a random 
matrix ensemble defined by V^(A) satisfy the orthogonal- 
ity relation 



dAe-^^(^)p„(A)P„(A) -<5„ 



(21) 



Then the exact expression for arbitrary N for the density 
of eigenvalues is 



N-l 



p(A)=e-^^(^)l^P„^(A) 



(22) 



i=0 
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All higher correlation functions can be recovered from 
the kernel 



i=0 



For the Cauchy random matrix ensemble the orthogo- 
nality relation ( pT| ) reads 



dA(A^ + l)-2A^P„(A)P,„(A) =5„, 



(24) 



We see that in contrast to the classical case only a fi- 
nite number of orthogonal polynomials exist. These are 
explicitly given by Jacobi polynomials analytically con- 
tinued to complex parameters |^. Indeed 



'{l+n-2N)n ^r(7V- 



22"n! 



r{N-n) 



{ix) 



(25) 



A second surprise comes from the fact that the eigenvalue 
distribution (E3) is exactly equal to 



7rA2 + l 



(26) 



and does not depend on N . There are no finite N correc- 
tions whatsoever to the spectral distribution. In particu- 
lar the classical short distance oscillations in the spectral 
density characteristic for RMM are absent. The Cauchy 
ensemble has, however, nontrivial iV-dependent 2-point 
correlation functions. 

The Levy-Smirnov ensemble can be analyzed starting 
from the CG distribution (111) with the measure (p^) 



1^ 



(27) 



A change of variables A; = 1 / Xi leads to 

Yldx,{e-^-')l[{x,^x,f (28) 



This is readily analyzed in terms of Laguerre polynomi- 
als similar to Chiral Unitary Ensembles (chGUE) 
Indeed the appropriate polynomials are 



Pn{x) = VNLUNx) 



(29) 



The eigenvalue density (^2|) can be rewritten, using the 
Christoffel-Darboux identity, as 



pix)^Ne-^^{L%_,iNx)L]^_,{Nx) - L%{Nx)L],_^{Nx)) 

(30) 

In particular we note that there exist a well defined mi- 
croscopic limit which corresponds to expressing the eigen- 
value density in terms of x = s/N"^ (i.e. on the scale of 
the eigenvalue spacing). The relevant eigenvalue density 
{l/N)p{s/N'^) is then given by 



p(s) = J2(2V^) + j2(2Vi) 



(31) 



Going back to the original variables the microscopic re- 
gion corresponds to the region of large eigenvalues A 
iV^A in the power-like tail. For these large eigenvalues 
we therefore observe chGUE-like oscillations 



P(A) 



1 

A2 
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J 



VA 



(32) 



Moreover we expect these oscillations to be universal in 
the following sense. A generic modification of the LS 
potential of the form 



nA)-21ogA+^ + || + ^ + 



El 

A3 



(33) 



will not change the oscillation pattern (^) . This follows 
from the results in [|l2| after the change of variables A — > 
X — 1/A. The coefficient in front of the logarithm cannot 
be changed, for otherwise the eigenvalue support becomes 
finite and the power like tails disappear altogether. The 
coefficients of the 1/A term (and higher gis) only affect 
the length scale of the universal oscillations. 

In the general case when the potential is of the form 
(pO|), and the asymptotic behavior of p{\) is like l/A^"*"", 
the mapping A — > a; = 1/A gives an effective potential 



V{x) - -2-Yieh 
a 



(34) 



General 



and an eigenvalue distribution p{x) ~ a;"' 
arguments [l^ show that the resulting eigenvalue spacing 
I/AtV" yields a microscopic distribution in the limit of 
iV — > oo with s = xN^^" fixed. The pertinent orthogonal 
polynomials should satisfy 



dxe 



(35) 



All this is very reminiscent of multicritical microscopic 
scaling and universality Jl3|-p^ in the classical random 
matrix case but here the "multicritical classes" are la- 
beled by a real parameter a and not by an integer. A 
thorough investigation of this new regime seems to be 
very interesting. 



V. SUMMARY 



^ After completing the paper, we noticed that a similar con- 
struction was recently used in |lO|. 



We have explicitly constructed matrix realizations of 
free random variables, with potential applications to a 
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number of stochastic phenomena. This opens several new 
venues for applying FRV calculus to Levy processes, in- 
cluding convolution, multiplication and addition of de- 
terministic matrix-like entries and other generalizations. 
Using the Coulomb gas analogy, we have shown that the 
exact matrix measure in the case of power-like spectra is 
non-local (involves determinants). The construction ex- 
hibits several non-trivial and novel features, among which 
the most interesting ones are a universal behavior in the 
tails of the distributions and an unusual large TV scaling. 
The expected microscopic eigenvalue distribution defines 
a new universal regime and represents a generalization 
of the multicritical scaling discussed in [p^-p5[ . Several 
of these issues, as well as practical applications of our 
results will be discussed in subsequent work. 
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